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Abstract 

We study arithmetic properties of short uniform random walks in arbitrary dimensions, 
with a focus on explicit (hypergeometric) evaluations of the moment functions and probability 
densities in the case of up to five steps. Somewhat to our surprise, we are able to provide 
complete extensions to arbitrary dimensions for most of the central results known in the two- 
dimensional case. 


1 Introduction 

An n-step uniform random walk in starts at the origin and consists of n independent steps 
of length 1, each of which is taken into a uniformly random direction. In other words, each step 
corresponds to a random vector uniformly distributed on the unit sphere. The study of such walks 
originates with Pearson [Pea05], who was interested in planar walks, which he looked at [Pea06] as 
migrations of, for instance, mosquitos moving a step after each breeding cycle. Random walks in 
three dimensions (known as random flights) seem first to have been studied in extenso by Rayleigh 
[Rayl9], and higher dimensions are touched upon in [Wat41, §13.48]. 

This paper is a companion to [BNSWll, BSW13] and [BSWZ12], which studied the analytic 
and number theoretic behaviour of short uniform random walks in the plane (five steps or less). In 
this work we revisit the same issues in higher dimensions. Somewhat to our surprise, we are able 
to provide complete extensions for most of the central results in the culminating paper [BSWZ12]. 

Throughout the paper, n and d are the positive integers corresponding to the number of steps 
and the dimension of the random walk we are considering. Moreover, we denote with v the half¬ 
integer 

^=^- 1 . ( 1 ) 

It turns out that most results are more naturally expressed in terms of this parameter and so we 
denote, for instance, with Pn{^‘i x) the probability density function of the distance to the origin after 
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n random unit steps in In Section 2.1, we mostly follow the account in [Hug95] and develop 
the basic Bessel integral representations for these densities beginning with Theorems 2.1 and 2.10, 
which are central to our analysis. In particular, a brief discussion of the (elementary) case of odd 
dimensions is included in Section 2.2. 

In Section 2.3, we turn to general results on the associated moment functions 

POO 

Wniiy;s)= / x''pn{v,x)dx. (2) 

Jo 

In particular, we derive in Theorem 2.18 a formula for the even moments 2k) as a multiple sum 

over the product of multinomial coefficients. As a consequence, we observe another interpretation of 
the Catalan numbers as the even moments of the distance after two random steps in four dimensions, 
and realize, more generally, in Example 2.22 the moments in four dimensions in terms of powers of 
the Narayana triangular matrix. We shall see that dimensions two and four are privileged in that 
all even moments are integral only in those two dimensions. 

In Section 3, we turn to a detailed analysis of the moments of short step walks: from two-step 
walks (§3.1) through five steps (§3.4). For instance, we show in (37), Theorem 3.6 and Example 3.25 
that the ordinary generating functions of the even moments for n = 2,3,4 can be expressed in terms 
of hypergeometric functions. We are also able to give closed forms for all odd moments for less 
than five steps. 

In Section 4, we perform a corresponding analysis of the densities Pn{^] x)' from two-step walks 
(§4.1) through five steps (§4.4). One especially striking result for n = 3, shown in Corollary 4.6, is 
the following functional equation for the probability density function p^{v]x). For 0 ^ x ^ 3, and 
each half-integer ^ ^ 0, the function F(x) := p'i(y\x)/x satisfies the functional equation 


F{x) 



3 — X 
1 -I- X 


(3) 


Finally, in Section 5, we make some concluding remarks and leave several open questions. 

As much as possible, we keep our notation consistent with that in [BNSWll, BSW13], and 
especially [BSWZ12], to which we refer for details of how to exploit the Mellin transform and similar 
matters. Random walks in higher dimensions are also briefly discussed in [Wanl3, Chapter 4]. In 
particular. Wan gives evaluations in arbitrary dimensions for the second moments, which we consider 
and generalize in Example 2.20, as well as for two-step walks. 


2 Basic results from probability 


2.1 The probability densities 

For the benefit of the reader, we briefly summarize the account given in [Hug95, Chapter 2.2] of how 
to determine the probability density function pn{v]x) of an n-step random walk in d dimensions. 
The reader interested in further details and corresponding results for more general random walks, 
for instance, with varying step sizes, is referred to [Hug95]. 

Throughout the paper, the normalized Bessel function of the first kind is 


jy{x) = v\ 



(-xV4)™ 
m\{m -\- 1 >)\ 


(4) 
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With this normalization, we have ji, (0) = 1 and 





/2 \ <^+ 1/2 
UJ 




(5) 


as a; —)■ c» on the real line. Note also that ji/ 2 ix) = sinc(x) = sin(x)/x, which in part explains why 
analysis in 3-space is so simple. More generally, all half-integer order ju{x) are elementary. 

Theorem 2.1. (Bessel integral for the densities, I) The probability density function of the 
distanee to the origin in d ^ 2 dimensions after n ^ 2 steps is, for x > 0, 

c\—i/ POO 

Pn{v;x) = — (6) 

Jo 

where, as introduced in (1), v = | — 1- 

Proof. Let JC be a random vector, which is uniformly distributed on the unit sphere in That 
is, X describes the displacement of a single step in our random walk. Then the Fourier transform 
of its induced probability measure p,x = ^ (II® II ~ 1) is 


iJx{q)= f = ji/(||q|i)- 


This is a special case of the famous formula [Hug95, (2.30)], 

/* /*oo / 9 \ ^/2~1 

y ^ e*i^''*>/(|lx|l)dx = 27r‘^/2y j Jd/2-i{tq)t‘^~^f{t)dt, (7) 

with q = ||q||, for integrals of orthogonally invariant functions. 

Note that the position Z of a random walk after n unit steps in is distributed like the sum 
of n independent copies of X. The Fourier transform of fj,z therefore is 


Az(g) = j”(llg||)- 


( 8 ) 


We are now able to obtain the probability density function pn{v,x) of the position after n unit 
steps in via the inversion relation 


= e fiz{q)dq. 

JR'I 

Combining (8), (9) and (7), we find 


(9) 


Pn{v,x) = 


(*00 


{21,) 


ly+i 


-Ju{tx)j{f{t)dt. 


Since the surface area of the unit sphere in is 27r'^+^/zz!, the density functions for the position 
and distance are related by 


Pn{v,x) = 


27r^+i||a;P‘"+i 


Vn{v\ Il®|| 


whence we arrive at the formula (6). 


( 10 ) 

□ 
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The probability densities psl^'x) of the distance to the origin after three random steps in 
dimensions 2, 3,..., 9 are depicted in Figure 1. With the exception of the planar case, which has 
a logarithmic singularity at x = 1, these functions are at least continuous in the interval [0,3], on 
which they are supported. Their precise regularity is provided by Corollary 2.5. For comparison, 
the probability densities of four-step walks in dimensions 2,3, ...,9 are plotted in Figure 2, and 
corresponding plots for five steps are provided by Figure 4. 




Figure 2: x) for ^ = 0, 5 ,1,..., | 

Observe how the density functions center and spike as the dimensions increase. Indeed, this is 
a general phenomenon and the distributions described by the densities Pn(y',x) approach a Dirac 
distribution centered at ^Jn. The intuition behind this observation is as follows: as the dimension 
d increases, the directions of each of the n random steps increasingly tend to be close to orthogonal 


4 














to each other. That is, given an incoming direction, the direction of the next step will probably 
belong to a hyperplane which is almost orthogonal to this direction. Pythagoras’ theorem therefore 
predicts that the distance after n steps is about ^/n. A precise asymptotic result, which confirms 
this intuition, is given in Example 2.4. 

Integrating the Bessel integral representation (6) for the probability density functions 
we obtain a corresponding Bessel integral representation for the cumulative distribution functions, 

P„(i/;x)=/ pn{^]y)dy, (11) 

Jo 

of the distance to the origin after n steps in d dimensions. 

Corollary 2.2. (Cumulative distribution) Suppose 2 and n ^ 2. Then, for x > 0, 


9 - 1 / foo 1 , 

P„(i/;x) = — / (tx)'"+Vj.+i(te)j”(t) —. 

Jo ^ 


( 12 ) 


Example 2.3. (Kluyver’s Theorem) A famous result of Kluyver [Klu06] asserts that, for n ^ 2, 


P„( 0 ; 1 ) 


1 

n + 1 


That is, after n unit steps in the plane, the probability to be within one unit of the starting point 
is l/(n + 1). This is nearly immediate from (12). An elementary proof of this remarkable result 
was given recently by Bernard! [Berl3]. 

It is natural to wonder whether there exists an extension of this result to higher dimensions. 
Clearly, these probabilities decrease as the dimension increases. 

(a) In the case of two steps, we have 




l+iy 




2 ^ r(z/ + 3/2) 
which, in the case of integers v 0 and x = 1, reduces to 

3fe 

3 47r^^(2/c+T)(f) 

Alternatively, 

In particular, in dimensions 4, 6 and 8, 


k-i/2nk? 

T{2k)' 


1 


1 3V3 




^2(3;!) = -- 


1 9^3 


3 207r 


(13) 


(14) 


and it is obvious from (13) that all the probabilities are of the form | for 

some rational factor Cj, > 0. 
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(b) In the case of three steps, we find 




p.o,lV 1 2048 

4 9457r2’ 


Indeed, for integers v ^ 0, we have the general formula 

^3(^; I) = z - ^ E - 3) 


/c=l 


r(2fc)r(3fc)' 


In the limit we arrive at the improbable evaluation 


5 F 4 


— 1 I I 1 

_§_ i 3 5 
11 ’ 3 ’ 2 ’ 3 


16 


= -TT 

16 


( 15 ) 


(16) 


since Psi^; 1) goes to zero as the dimension goes to infinity. 

The case of n-step walks, with n > 4, is much less accessible [BBBG08, §5], and it would be 
interesting to obtain a more complete extension of Kluyver’s result to higher dimensions. o 

Example 2.4. Asymptotically, for x > 0 as —>■ 00 , 


ju{t) ~ exp 


—)■ 

4z/ + 2 y 


It may thus be derived from the Bessel integral representation (6) that asymptotically, as the 
dimension goes to infinity, Pn{v;x) ~ qn{v;x), where 


qn{i^\x) 


2 "^' /2i^+l 
vl \ n 


I'+i 


exp 


2 i/+ 1 

2n 


for X ^ 0 and g„(j^;x) = 0 for x < 0. The probability density qn{v]x) describes a scaled chi 
distribution with 2^ + 2 = d degrees of freedom. Its average is 



xqniy] x)dx 


2ti ^ 4” 

2j^ + 1 r(z/ +1) ’ 


which converges to y/n as —>■ 00 . More generally, for the sth moment, with s > — 1, 



x^qn{v, x)dx 


2n Y^^r{u+^ + l) 

r(^^ + i) 


s(s + 2 ) 
8 z^ 


2z/+ 1 
2nv 
2v+l 


s/2 


1 + 




as V ^ 00 , and it is straightforward to compute further terms of this asymptotic expansion. The 
fact that the sth moment approaches for large dimensions, of course, reflects the observation 
from Figures 1, 2 and 4 that the probability densities Pn{v\x) center and spike at approximately 
i/n. o 
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These Bessel integral representations also allow us to deduce the regularity of the density func¬ 
tions. 


Corollary 2.5. (Regularity of the density) The density pn(v;x) is m times continuously 
differentiable for x > 0 if 


n > 


m 1 

V 1/2 


+ 1 ) 


or, equivalently, m < in — l){v 1/2) — 


Proof. It follows from ( 6 ), that 


Pnit^^x) 


n —2i/ roo 


(17) 


Observe that this integral converges absolutely il 21/ + 1 — {n + l){v + 1/2) < —1, in which case 
Pn{v',x) is continuous. Repeatedly differentiating under the integral sign as long as is permitted, 
we conclude that pn(y]x) is m times continuously differentiable if 


2j^ -I- 1 — (n -b l)(i^ -I- 1/2) < —1 — TO, 


and it only remains to solve for n, respectively to. □ 

Example 2.6. In the case d = 2, this implies that pn{0',x) is in for n > 3, in for n > 5, in 
for n > 7, and so on. Indeed, p3{0-,x) has a logarithmic singularity at x = 1, and p4{0;x) as 
well as P5(0; x) are not differentiable at x = 2,4 and x = 1, 3, 5, respectively. o 

Corollary 2.7. Let n ^ 4 and d^ 2 such that (n, d) ^ (4,2). Then, 

where the derivative is understood to be taken from the right. 

Proof. Starting with (17), we differentiate 2v + 1 many times, and compare with 


Pn-i{v;l) 


0 — 2 h' 


t 


2i^+l An 


j”(t)dt. 


□ 


In the case d = 2, this reduces to p(j(0;0) = p„_i(0;l), which was crucial in [BSWZ12] for 
explicitly evaluating ^ 5 ( 0 ; 0 ). 

Example 2.8. As long as p„(t^;x) is sufficiently differentiable at x = 1, we can further relate the 
value p„(i/;l), occuring in (18), to corresponding values of derivatives. Let us illustrate this by 
showing that 

, , , 2nv + n — 1 , . , , 

Pn{v\ 1) = -- Pn{v-, 1 ) (19) 

for all n ^ 3, i/ > 0 such that (n, v) 7 ^ (3,1/2). In that case, we may, as in Corollary 2.5, differentiate 
(17) under the integral sign to obtain 

o —2I2 poo 

Pni^l 1) = (2i^ + l)p„(!^; 1) + t^''^^jl{t)j//{t)dt 
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On the other hand, integrating (17) by parts, we find 


1 ) = - 


n + 1 2 
21 ^ + 2 ^!^ 


— 2 i/ /-oo 




Combining these, we arrive at (19). o 

The densities of an n-step walk can be related to the densities of an (n — l)-step walk by the 
following generalization of [Hug95, (2.70)]. 

Theorem 2.9. (Recursion for the density) For x > 0 and n = 1, 2,..., the function 

v\ Pn{v;x) 


ifniy^x) = 


2'k’'+^ 


satisfies 


^\2n2L' rl 

lfn{v]x) = ' / fin-liv, \/l + 2Xx + X‘^){1 - 

[Zvy.TT 


( 20 ) 


IS 


Proof. Recall from (10) that the probability density of the position after n steps in 

Pn{v,x) = ifniv] Ila:|l). 

Since the steps of our walks are uniformly distributed vectors on the unit sphere in we have 

Pn(y\x)= / pn-i(v]x-y)dS = / ij}n-i{vA\x-y\\)dS, 

where dS' denotes the normalized surface measure of the unit sphere. After introducing d-dimensional 
spherical polar coordinates, as detailed, for instance, in [Hug95, p. 61], and a straightforward change 
of variables we arrive at (20). □ 


Finally, a computationally more accessible expression for the densities Pn{h';x) is given by the 
following generalization of a formula, which was derived by Broadhurst [Bro09] in the case of two 
dimensions, that is, z/ = 0. Note that (6) is the special case fc = 0 in (21). 

Theorem 2.10. (Bessel integral for the densities, II) Let n ^ 2 and d 2. For any 

nonnegative integer k, and x > 0, 

n-v 1 fOO / 1 J \ * 

Pn(j^;a:) = {txfi+'^+^J^+kitx) j”(t)dt. (21) 

Proof. As in [Bro09], we proceed by induction on k. The case A: = 0 is (6). Suppose that (21) is 
known for some k. By the Bessel function identity 

^{z''Mz))=z''J,.i{z) ( 22 ) 

az 

we find that, for any smooth function g{t) which is sufficiently small at 0 and oo, 


{tx)°‘~'~^Ja{tx)g{t)dt = — 


— (te)“+Va(te) 
at 


g{t)dt 


1 r°° 

= — / (te)“+V„(te) 

X Jo 

1 poo 

= — (te)“+^Ja(te) 

Jo 


d , . 


dt 


Id 


dt. 












In the second step, we used integration by parts and assumed that g(t) is such that 

{tx)°‘~^^Jaitx)g{t) (23) 

vanishes as t —>■ 0 and t —>■ oo. In the present case, a = v + k and g{t) = (—j"(0- order to 
complete the proof of (21), it only remains to demonstrate that (23) indeed vanishes as required. 

As t —?> 0, we have {tx)°‘'^^ Ja{tx) = C)(t^“+^) = while g{t) = 0(1), because j^{t) is 

an even function. Hence, for i/ > —1, the term (23) indeed vanishes as t —)■ 0. On the other hand, 
as t —)■ oo, we have (ta;)“+^ Ja(te) = 0(t“+^/^) = 0(t'^+^+^/^) by (5). Moreover, (5) and (22) imply 
that = 0(t“'^“^/^), as t —)■ oo, and therefore that g(t) = 0(t“"(‘"+i/2)-fc). Since n > 1, it 

follows that (23) also vanishes as t —>■ oo. □ 

We note that the (—yay)^j"(i) in the integrand of (21) may be expressed as a finite sum of 
products of Bessel functions. This is made explicit in Remark 2.19. 


2.2 The probability densities in odd dimensions 

In the case of odd dimension d, the Bessel functions in Section 2 have half-integer index i/ and are 
therefore elementary, so that the situation is fairly well understood since Rayleigh [Rayl9, Hug95]. 
In particular, the probability density functions Pn{v] x) are piecewise polynomial in odd dimensions. 
This is made explicit by the following theorem, which is obtained in [GP12], translated to our 
notation. 


Theorem 2.11. (Density in odd dimensions [GP12, Theorem 2.6]) Assume that the di¬ 
mension d = 2m 1 is an odd number. Then, for 0 < x < n, 


p„(m - ^;x) = 


(2a;)2™r(m) 

r(2m) 


^_d 

2x dx 


P, 


i(x), 


where Pm,nix) is the piecewise polynomial obtained from convoluting 


r(m-|-l/2) f (1 — a;^)™ 1, for X G [—1,1], 
r(l/2)r(m) ( 0, otherwise, 


n — 1 times with itself. 

Example 2.12. In the case n = 3 and d = 3, we have to = 1 and 

P 1 / 1. if 

2 10, otherwise, 


as well as 


and, hence, 


p t \ p (AP / PAe ^ f 2-|a;|, if |a;| < 2, 

Pi. 2 (x) = J_ P,iim,dx - t)dt = - I otherwise, 


3 — x^, 

P,3ix) = I Pi, 2 it)Pi,i{x - t)dt = ^ ^ i(|a:|-3)2, ifl<|a:|^3, 

0 , 


if \x\ ^ 1, 
if 1 < \x\ s 
otherwise. 
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We thus find that, for 0 < x < 3, 


P3{\]x) = 


---) 
2x dx j 


-Pi.3 (a;) 


X f 2x, if 0 ^ a; < 1 , 
4 \ 3 — X, if 1 < X < 3. 


Since ji/ 2 (x) = sinc(x) = sin(x)/x, evaluation of the densities in three dimensions can also be 
approached using the tools provided by [BBOl]. 

Similarly, we obtain, for instance. 


,i X _ a; J x (8 —3x), if 0 < x < 2 , 
“ 16 \ (4 - x)2, if 2 < X ^ 4. 


(24) 


We note that Theorem 2.11 can be usefully implemented in a computer algebra system such as 
Maple or Mathematica. o 

Example 2.13. (Moments in odd dimensions) By integrating (24), we are able to symbolically 
compute the corresponding moments, as introduced in ( 2 ), as 


Wi{l/2-s) 


= 2 '*+^ 


1 _ 2'>+2 1 
(s + 2) (s + 4) (s + 3) 


which has a removable singularity —2 and poles at —3 and —4. Likewise, in five dimensions, 

(12)32"+! (s3 + 27 s2 + 230 s + 616 + 64 2« (s3 + 15 + 62 s + 56)) 

^ ’ (s + 12) (s + 10) (s + 9) (s + 8) (s + 7) (s + 6 ) (s + 5) (s + 4) (s + 2) 


with poles at —5,—7,— 8 ,—9,—10,—12 and removable singularities at the other apparent poles. 
Thus, even this elementary evaluation has subtle structure. o 


2.3 The moment functions 


Theorem 2.10 allows us to prove a corresponding Bessel integral representation for the moment 
function 


Wn{i';s)= / x®p„(z/;x)dx 

Jo 


of the distance to the origin after n random steps. The following result generalizes [Bro09, Theo¬ 
rem 1 ] from two to arbitrary dimensions. 

Theorem 2.14. (Bessel integral for the moments) Let n > 2 and d ^ 2. For any nonnegative 
integer k, 


Wniv; s) = 


2,-fc+ir (I 1) 


r(z. + i)r(fc-f) Jo 

provided that k — n{i' -I- 1 / 2 ) < s < 2k. 

Proof. Using Theorem 2.10, we have 


X dx / 


( 25 ) 


c\ — i/ poo p 

v'- Jo Jo 


(txY+’^+^J^+kitx) ( ) j”(t)dMx. 
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Interchanging the order of integration and substituting z = ta;, we obtain 


Wn{iy;s) = 


±2k—s — l 


1 d 
t dt 


m 


^v+s — k+l 


Jj.+fc(z)dzdt. 


The inner integral may be evaluated using the standard Bessel integral evaluation 


2 “r ( 

z“J^(z)dz = -jrA 


ap f l+a+i>' ^ 


( 1 ^) 


which holds for a and z/ such that a + > — 1 and a < 1/2. We conclude that 


Wniv;s) = 


-y—y poo 


±2k-s-l 


-is' 


2‘^+"-'=+ir(f + 1/ + 1) 
r (aT- f) 


dt, 


which is the desired result. For the evaluation of the intermediate Bessel integral, we assumed 
s > —2v — 2 = —d and s < k — v — 1/2, and so (25) holds for all s in this non-empty strip provided 
that the original integral converges. Using the asymptotic bounds from the proof of Theorem 2.10, 
we note that the integral (25) converges absolutely in the strip k — n{v + 1 / 2 ) < s < 2k. Analytic 
continuation therefore implies that (25) holds for all s in this strip. □ 

We deduce the following from (25), with k = 0. This extends [BSW13, Proposition 2.4]. In 
particular, we observe that, for n> 2, the first pole of Wn{i^] s) occurs at s = —{2i' -I- 2 ) = —d. 

Corollary 2.15. (Poles and residues of the moments) Let n > 2. In the half-plane Re s > 
—n{v 1/2), the moment functions Wn{v', s) are analytic apart from simple poles at s = —d — 2m 
for integers m such that — i) — The residues of these poles are 


R 0 S 5 —2m ^^n(^7 '5) — 


2^ — 2u—2m ^ 

m)\ ml 


^2v-\-2m+l 


(a;)da 


Proof. By equation (25), with fc = 0, we have 


Wn{v;s) 


2'*+ir(f i/-f 1) 

r(z. + i)r (-§) 


'j”(^)da 


valid for s in the strip —n{v -|- 1/2) < s < 0, in which the integral converges absolutely. In the 
region of interest, the only poles are contributed by the factor T (| -|- z/ -|- l), which, as a function 
in s, has simple poles at s = —d — 2m, for m = 0 , 1 , 2 ,. .., with residue 2 *' . □ 

Note that the value for the residue of lU„(z^; s) at s = —d agrees with 


from (18) provided that these values are finite. 
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Example 2.16. The moment functions W 3 {i';s) are plotted in Figure 3 for 1 ^ = 0, 1,2. The first 
pole of W 3 {v\ s) occurs at s = —d, is simple and has residue 

2”2i/ roo 2 

Res^=-dW3(i^;s) = / x'^''+'^jl{x)dx = 

do VdTT j 

as follows from the Bessel integral formula (72). In order to record some more general properties 
of the pole structure of W'i{v\ s), we use the fact that W3(^; s) satisfies a functional equation, (48), 
which relates the three terms lT3(z^; s), W^{v] s + 2) and Wz{v] s + 4). By reversing this functional 
equation, we find that the residues = ReSs=_d-2fe W^(v] s) satisfy the recursion 


9(fc + l)(fc -\- V -\- l)rfc_|_i 
\ ^20 {k + — 20 (fc + i) + 1^ rfe — (fc — v)(k — 2i/)ri. 


-1) 


(26) 


with r_i = 0 and xq = Observe that, when i/ = 0, the recursion for these residues is 

essentially the same as the recurrence (48) for the corresponding even moments (with Uk replaced 
by 'd^^Tk)- As recorded in [BSW13, Proposition 2.4], this lead to 


Ress=_2(/c+i) W3(0;s) = 


2 lF3(0;2fc) 

TStt 


Define, likewise, the numbers V 3 {v]k) by 


2 3" V3{,^-,k) 

dieSs=-2{v+k + l)W^{v,s) — 

In analogy with (26), we find that Uk = k) solves the three-term recurrence 
(fc -|- l)(fc -\- V -\- l)Ufc+l 

= \ ^20 (k + — 20(k+\)v — 4i/^ -|- 1^ mj, — 9(fc — v){k — 2v)uk-i, 

with M_i = 0 and uq = 1. For small dimensions, initial values for k) are given by 


d = 2 (i/ = 0) 

d = 4 iiy=l) 
d = 6 {v = 2) 
d = 8 (v = 8) 


1,3,15,93,639,4653,35169,272835,2157759,... 

1, -2, -2, -6, -24, -114, -606, -3486, -21258,... 
1,-5,6,2,6,18,66,278,1296,... 

1, -f, 21, -20,0, -9, -20, -60, -210,... 


Note the increasingly irregular behaviour of V^iy^k) as d increases. In particular, in dimension 8, 
we find V3(3; 4) = 0, which signifies the disappearance of the perhaps expected pole of 114(3; s) at 
s = —16. o 


Example 2.17. In the case n = 4 and = 1, that is d = 4, the moment function 114(1; s) has a 
simple pole at s = —4 with residue 


ReSs=_4 114(1; s) 


1 

4 



x^ji{x)dx = 


4 
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Figure 3: W^{v\ s) on [—9, 2] for i/ = 0,1, 2. 


and is otherwise analytic in the half-plane Re s > —6. At s = —6, on the other hand, W^ll; s) has 
a double pole. Indeed, analyzing the functional equations that arise from (63) and (65), we derive 
that 

lim (s -I- 6)^W4(1; s) = lim (s + 2)^W4(0; s) = -7^, 

where in the last equality we used the known value from the planar case [BSWZ12, Example 4.3]. 
Similarly, we obtain 

ReSs=_6 W4(l;s) = -7 ReSs=_2 W4(0; s)-h-*7 lim (s-h 2)^W4(0; s) 

2 24 s—2 

1 91og(2) 

167r^ 47r^ 

In the higher-dimensional case, W 4 (iy; s) has poles at s = —d — 2m for m = 0,1, 2,..., which are 
initially simple but turn into poles of order (up to) 2 beginning at s = —{Ai/ -1-2). o 

The approach indicated in Examples 2.16 and 2.17 enables us to determine, at least in principle, 
the pole structure of the moment functions s) in each case. We do not pursue such a more 

detailed analysis herein. 

We next follow the approach of [Bro09] to obtain a summatory expression for the even moments 
from Theorem 2.10. 


Theorem 2.18. (Multinomial sum for the moments) The even moments of an n-step random 
walk in dimension d are given by 


Wn{v,2k) = 


(k + v)\v\'^-^ 
(k -\- nv)\ 


E 


ki-\ - \-kn—k 

Proof. Replacing fc by fc -I- 1 in (25) and setting s = 2k, we obtain 


ki,... ,knj \ki + 1 ^,... ,kn + 


Wn{m,2k) = 




ik + u)\ f 2d 


-;a5 1 


21=0 


Observe that, at the level of formal power series, we have 


) E 

X da: ' 




m^O 


= klok- 


x—0 


(27) 
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Recall from (4) the series 


to conclude that 


j^(x) 


m^O 


(-^2/4)m 

m\{m + v)\' 


Wnii^;2k) 


{k + i/)! 


v\^k\ 


E 


miH- \-mn—k 


mi! • • • m„! (mi + ly)! ■ ■ ■ (rUn + J^)! ’ 


which is equivalent to the claimed formula. 


□ 


Remark 2.19. Proceeding as in the proof of Theorem 2.18, we observe that 



(y + k)\ 


jv+k (3 


and hence 



k 

>1 (a;) • • ■ (x) 


= E 

fcl-|- \-kn — k 


kl _ i^i! ■ ■ • _ 

fell • • • fc„! (fci + i/i)! ■■ ■{kn + jynV- 




If applied to (27), this (finite) expansion, together with ju{0) = 1, makes the conclusion of The¬ 
orem 2.18 apparent. In conjunction with the asymptotics of ji,, we conclude, as in the proof of 
Theorem 2.10, that the integrand in (21) is as t —>■ 00 , so that each additional 

derivative improves the order at -l-oo by 1 — at the expense of increasing the size of the coefficients. 
We note, inter alia, that (27) may, alternatively, be expressed as 


k^O 


Wnjiy; 2k) 
k\{k + v)\ 



which yields a fine alternative generating function for the even moments. 

Example 2.20. In the case fc = 1, Theorem 2.18 immediately implies that the second moment of 
an n-step random walk in any dimension is 


W'„(i/; 2) = n. 


This was proved in [WanlS, Theorem 4.2] using a multi-dimensional integral representation and 
hyper-spherical coordinates. Similarly, we find that 


Wniiy-A) 


n{n{v -b 2) — 1) 
v+l ' 


(28) 


More generally. Theorem 2.18 shows that Wn{v\ 2k) is a polynomial of degree k in n, with coefficients 
that are rational functions in v. For instance, 


r ri{n^{iy + 2){iy + i)-3n{u + 3) + A) 

(Z.+ 1)2 


(29) 


and so on. 


o 
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Using the explicit expression of the even moments of an n-step random walk in dimension d, we 
derive the following convolution relation. 

Corollary 2.21. (Moment recursion) For positive integers ni,n 2 , half-integer ly and nonneg¬ 
ative integer k we have 

W^ni+n2(^^; 2fc) = ^ J - j))- (30) 

Note the special case n 2 = 1, that is 


W„{iy;2k)=Y, 

j=o 


[k + vy-vl 


jj (fc - j + 


H"n-i(u2j), 


(31) 


which allows us to relate the moments of an n-step walk to the moments of an (n — l)-step walk. 

Example 2.22. (Integrality of two and four dimensional even moments) Corollary 2.21 
provides an efficient way to compute even moments of random walks in any dimension. For illus¬ 
tration, and because they are integral, we record the moments of an n-step walk in two and four 
dimensions for n = 2, 3,..., 6. 


lU 2 ( 0 ; 2 fc) 

fU 3 ( 0 ; 2 fc) 

W4(0;2fc) 

lU 5 ( 0 ; 2 fc) 

lU 6 ( 0 ; 2 fc) 


1 . 2 . 6 . 20 . 70 . 252 . 924 . 3432 . 12870 .. .. 

1 . 3 . 15 . 93 . 639 . 4653 . 35169 . 272835 . 2157759 .. .. 

1 . 4 . 28 . 256 . 2716 . 31504 . 387136 . 4951552 . 65218204 .. .. 

1 . 5 . 45 . 545 . 7885 . 127905 . 2241225 . 41467725 . 798562125 .. .. 

1 . 6 . 66 . 996 . 18306 . 384156 . 8848236 . 218040696 . 5651108226 .. .. 


For n = 2, these are central binomial coefficients, see (36), while, for n = 3,4, these are Apery-like 
sequences, see (40) and (61). Likewise, the initial even moments for four dimensions are as follows. 


fU 2 (l; 2 fc) 
W^ 3 (l; 2 fc) 
1 ^ 4 ( 1 ; 2 fc) 
W^{l-2k) 
W^{\-2k) 


1, 2, 5,14,42,132,429,1430,4862,... 

1.3.12.57.303.1743.10629.67791.448023.. .. 

1.4.22.148.1144.9784.90346.885868.9115276.. .. 

1.5.35.305.3105.35505.444225.5970725.85068365.. .. 

1.6.51.546.6906.99156.1573011.27045906.496875786.. .. 


Observe that the first terms are as determined in Example 2.20. In the two-step case in four 
dimensions, we find that the even moments are the Catalan numbers Cfe, that is 


W2{l-,2k) 


(2fc + 2)! 1 f2k\ 

(fc + l)!(fc + 2)! k + l\k)' 


(32) 


This adds another interpretation to the impressive array of quantities that are given by the Catalan 
numbers. 
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It is a special property of the random walks in two and four dimensions that all even mo¬ 
ments are positive integers (compare, for instance, (28)). This is obvious for two dimensions from 
Theorem 2.18 which, in fact, demonstrates that the moments 


W"n(0;2fc) 


E 

All H- \~^n —A; 




2 


count abelian squares [RS09]. On the other hand, to show that the four-dimensional moments 
Wn{^] 2fc) are always integral, it suffices to recursively apply (31) and to note that the factors 


'k\ {k + l)\ _ 1 fk\ /k +1 \ 

j) {k-j + l)!(j + 1)! j + 1 U7 Vi/ 


(33) 


are integers for all nonnegative integers j and k. The numbers (33) are known as Narayana numbers 
and occur in various counting problems; see, for instance, [Sta99, Problem 6.36]. o 

Example 2.23. (Narayana numbers) The recursion (31) for the moments Wn{v]2k) is equiva¬ 
lent to the following interpretation of the moments as row sums of the nth power of certain triangular 
matrices. Indeed, for given z/, let A(n) be the infinite lower triangular matrix with entries 




{kv)\v\ 

U7 {k - j + y)\[j + v)\ 


(34) 


for row indices fc = 0,1,2,... and column indices j = 0,1,2,.... Then the row sums of are 

given by the moments Wn+i{v-, 2k), k — 0,1,2,.... For instance, in the case n = 1, 


kl(l) 


'1 

0 

0 

0 



■ 1 

0 

0 

0 


1 

1 

0 

0 



3 

1 

0 

0 


1 

3 

1 

0 


, 7(1)3 ^ 

12 

9 

1 

0 


1 

6 

6 

1 


57 

72 

18 

1 














with the row sums 1, 2, 5,14,... and 1,4, 22,148,... corresponding to the moments 11^2(1; 2fc) and 
117(1; 2fc) as given in Example 2.22. Observe that, since the first column of A{v) is composed of 
I’s, the sequence of moments 117(1; 2k) can also be directly read off from the first column of 7(z/)". 
The matrix 7(1) is known as the Narayana triangle or the Catalan triangle [Slol5, A001263]. o 


Remark 2.24. Let us note another point of view on the appearance of the Narayana triangle in 
the context of random walks. Let X be a random vector, which is uniformly distributed on the 
unit sphere in If 9 is the angle between X and a fixed axis, then A = cosd has the probability 
density [Kin63, (2)] 


^{v-\/2)\ 


(1-a2)"-i/7 


A G [—1,1]. 


Denote with the random variable describing the distance to the origin after n unit steps. Then 
Rn+i is related to i?„ via [Kin63, (9)] 


Rn+\ — \/l -1 2Ai?„ -|- R'^. 
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Writing E[X] for the expected value of a random variable X, we therefore have 


2k) = E[Rl%,] = E[(l + 2Ai?„ + RI)% 


In terms of the generalized Narayana polynomials 

= E[(l + 2Ay^ + z)'=-i], (35) 

which were introduced in [AMV13, (6.2)] for A: ^ 1, we obtain 

W^+,im,2k) =E[Mt''^,iRl)]. 


This recurrence identity on the moments can then be expressed in matrix form by defining a matrix 
A{v) as in Example 2.23. 

We note, moreover, that expressing the Narayana polynomials (35) in terms of the Gegenbauer 
polynomials C'^{z), as demonstrated in [AMV13, Theorem 6.3], we deduce the expression of the 
even moments as 


Wn+i{v,2k) = 


k\ 


{2v + l)j. 


E 


n n2 / 1 + ^ 


which is a variation of (31). 


Example 2.25. (Six dimensional even moments) To contrast with the integral even moments 
in 2 and 4 dimensions in Example 2.22, we record a few initial even moments in 6 dimensions. 


1^2(2; 2fc) : 1,2,14/3,12,33,286/3,286,884,8398/3,... 

1^3(2; 2fc) : 1,3,11,139/3,216,1088,5825,32763,191935,... 

1T4(2; 2k) : 1,4,20, 352/3, 2330/3,16952/3,133084/3, 370752, 3265208,... 

It may be concluded from (34), with i/ = 2, that the entries of the matrix A{2) satisfy Ai, j{2) G 
This implies that the even moments in dimension 6 are rational numbers whose denominators are 
powers of 3. Similar observations apply in all dimensions but we do not pursue this theme further 
here. o 


3 Moments of short walks 


3.1 Moments of 2-step walks 


It follows from Theorem 2.18 that the general expression of the even moments for a 2-step walk in 
d dimensions is given by 


lT2(zz;2fc) 


cm 



u\{2k -\- 2u)\ 

{k + v)\{k + 2v)\' 


(36) 


Example 3.1. Note that in the special case of = 0, that is, dimension 2, this clearly reduces to 
the central binomial coefficient. In dimension 4, as noted in Example 2.22, the even moments are 
the Catalan numbers. We note that the generating function for the two-step even moments is 


Y,W2(,v-2k)x^ 


2 F 1 


2v + l 



(37) 
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which reduces to the known generating functions for v = 0,1. Equation (37) is an immediate 
consequence of rewriting (36) as W'2(^;2fc) = 2^*^ ^ 

In fact, (36) also holds true when k takes complex values. This was proved in [Wanl3, The¬ 
orem 4.3] using a multi-dimensional integral representation and hyper-spherical coordinates. We 
offer an alternative proof based on Theorem 2.14. 


Theorem 3.2. For all complex s and half-integer v ^ 0, 

iy\r{s-^2n + l) 


W2{iy;s) = 


r (f -I-1/ -I-1) r (I -I- 2i/ -I-1) 
Proof. The case fc = 0 of (25) in Theorem 2.14 gives 


w (...)- + - + Tx— S-(x)dx 


r(^. + i)r (-§) 

provided that —n{v -\- 1/2) — 1 < s < 0. Using that, for 0 < Re s < 2z/ -|- 1, 



^/^(a;)da: = 2^ 


r(f)r(df^ + ^)r(i + ^)^ 
r(i)r(i-f + ^)r(i-§ + 2i.)’ 


(38) 


the claimed formula then follows from the duplication formula for the gamma function and analytic 
continuation. □ 


3.2 Moments of 3-step walks 

Lemma 3.3. The nonnegative even moments for a 3-step walk in d dimensions are 


W:i{m,2k)=Y, 

3=0 


k -\- A f2j 2v\ [j -\- V 


J 


J 


J 


-2 


(39) 


Proof. We apply Corollary 2.21 with ni = 2 and 712 = 1- Using that lUi(^;2fc) = 1 and that an 
evaluation of 1^2(1^; 2k) is given by (36), we obtain 


WsA; 2k) 



(k -\- v)\v\ 

{k- j -\-n)\{j -\-i^)\ 


W2{m,2j) 


[k-\-u)\i3\ i3\[2j2v)\ 

[k- j + v)\{j + v)\ {j -\- ty)\{j -\- 2n)l' 


Expressing the factorials as binomial coefficients yields (39). □ 

Example 3.4. (Generating function for 3 steps in 2 dimensions) In the case d = 2, or 
13 = 0, the moments of a 3-step walk reduce to the Apery-like numbers 


1^3(0; 2fc) 



(40) 
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[BSWZ12, (3.2) & (3.4)] show that the generating function for this sequence is 


Y,W3iO-,2k)x'^ 

k=0 


1 


1 + 3a; 


2.Fi 


1 2 
3’ 3 
1 


27a;(l — a;)^\ 

(1 + 3a;)3 J ■ 


o 


Example 3.5. (Generating function for 3 steps in 4 dimensions) In the case d = 4, 
or = 1, the moments, whose initial values are recorded in Example 2.22, are sequence [Slol5, 
A103370]. The OEIS also records a hypergeometric form of the generating function (as the linear 
combination of a hypergeometric function and its derivative), added by Mark van Hoeij. On using 
linear transformations of hypergeometric functions, we have more simply that 


J^W3(l;2k)x^ 


1 1 p (\,i 21x{l-xf \ 

2 x^ X 2a;2(l + 3a:)^ 2 (l + 3a:)3 )' 


which we are able to generalize. o 

Example 3.5 suggests that a nice formula for the generating function for the moments W3(jz; 2k) 
exists for all even dimensions. Indeed, we have the following result. 

Theorem 3.6. (Ordinary generating function for even moments with three steps) For 

integers v ^ 0 and |a;| < 1/9, we have 


Y,W3{v;2k)x’^ 


(_l)^-(l_ 1/^)2. / 1,2 27 x { 1 - x )^ \ (1 

l + 3x ^ + (l + 3x)^ J ^ 


(41) 


where qv{x) is a polynomial (that is, qi,{l/x) is the principal part of the hypergeometric term on the 
right-hand side). 

Proof. For integers n ^ 0, define the rational numbers H{v, k) by 


(-ir ( |,i 21x{l-x)^ \ 

(2/') l + 3x ^ ^^1 + ^ (l + 3a;)3 ) 


^ H{v,k)x'" 

k— — 2v 


(42) 


for k ^ —2v, and H{v, fc) = 0 for fc < —2v. Writing the sum (39) for 2k) in hypergeometric 

form, we obtain, for integers k,v ^ 0, the representation 


W3{v,2k) 


3 F 2 


f —k,—k — v,v -\- 1/2 
\ + 1, 2v + 1 



(43) 


In order to prove the claimed generating function (41) it therefore suffices to show the (more precise) 
claim 


H{v, k) = Re 3E2 


f —k,—k — v,v -\- 1/2 
\ jz + 1, 2v + 1 



(44) 


For instance, this predicts that q 2 {x) = 1/6 — 5/6a; + + 1/3+ x'^. 

Taking the real part is only necessary for k < —v while, for k ^ —v, the 3F2 is terminating. 
We note, as will be demonstrated later in the proof, that the right-hand side of (44) vanishes for 
k < — 21 /, that is, the 3F2 takes purely imaginary values then. 
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The holonomic systems approach [Zei90], implemented in the Mathematica package Holonomic- 
Functions, which accompanies Koutschan’s thesis [Kou09], shows that the coefficients H{v,k), 
defined by (42), satisfy the recursive relation 

9(fc + l)(fc + V k) = \ ^20 {k + 1)^ + 60 (fc + I) + 1^ H{v,k + 1) 

-{k + 2v + 2){k + Zv + 2)H{v,k + 2) (45) 


for all integers k and all integers z/ ^ 0. We already know that (44) holds for 1 ^ = 0 and fc ^ 0. 
Verifying, by using 


3F2 


/ 1 , 1 , 1/2 
V 1.1 



1 


and letting x —> 1 to see that the real part vanishes (for any choice of analytic continuation to 
X = 1), that (44) holds for v = Q and k = —1, we conclude from (45) that (44) is indeed true for 
V = 0 and all integers k. 

As in the case of (45), we find that the coefficients H{v]k) further satisfy the dimensional 
relations 


2{k + 2u){k + 2,1' — l)(fc + 2>v)H{v, k) 

= v^{lk+ lbv — A)H{v—l]k+ 1)— ^v^{k+ v)H{v—l]k) (46) 

and 

2{k + l)(fc + 2v)H(v- k) = v^H{v - 1; fc + 2) - 2v^H{v - 1; fc + 1) (47) 

for all integers k and all integers v ^ The three relations (45), (46), (47), together with the 
case = 0 as boundary values, completely determine the coefficients H{u] k) for all integers k and 
z^ ^ 1. 

Since we already verified the case z/ = 0, it only remains to demonstrate that the right-hand side 
of (44) satisfies the same recursive relations. Another application of HolonomicFunctions finds 
that the 3 F 2 on right-hand side of (44), and hence its real part, indeed satisfy (45), (46), (47) for 
the required (real) values of i' and k. □ 

For the convenience of the reader, and because we will frequently use it in the following, we 
state Carlson’s Theorem next [Tit39, 5.81]. Recall that a function f[z) is of exponential type in a 
region if |/(z)| ^ for some constants M and c. 

Theorem 3.7. (Carlson’s Theorem) Let f be analytic in the right half-plane Re z ^ 0 and of 
exponential type with the additional requirement that 

\f{z)\^Me^'^' 

for some d < n on the imaginary axis Re z = 0. If f{k) = 0 for k = 0,1, 2,..., then f(z) = 0 
identically. 

Example 3.8. Applying creative telescoping to the binomial sum (39), we derive that the moments 
IT 3 (z^; 2k) satisfy the recursion 

{k + 2v + l){k + 2 ,v+l)W 3 {v] 2 k + 2) (48) 

= 4 (20 (fc -F \f + m{k+\)v + 36z/2 -f 1) 2k) - 9fc(fc -F v)W 3 {v] 2k - 2). 
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Observe that Wn{y',s) is analytic for Re s ^ 0 and is bounded in that half-plane by |W„(i^;s)| < 
nRe s (^because, in any dimension, the distance after n random steps is bounded by n). It therefore 
follows from Carlson’s Theorem, as detailed in [BNSWll, Theorem 4], that the recursion (48) 
remains valid for complex values of fc. o 


The next result expresses the complex moments W 3 {iy;s) in terms of a Meijer G-function and 
extends [BSW13, Theorem 2.7]. 


Theorem 3.9. (Meijer G form of W 3 ) 


For all complex s and dimensions d ^ 2, 


1 ^ 3 ( 1 /; s) = 


r(f -f 3,3 


1,1 + z/, 1 + 2iy 


2 ’ 



Proof. If 0 < Re ( 5 ) < i/ + |, then 


jij{x)dx = 2 


^_, r(f)r(i + ^) 
r(i-f + ^)- 


This simple integral is a consequence of the fact that Wi{v]s) = 1 combined with Theorem 2.14 
with n = 1 and fc = 0. Similarly, for 0 < Re (s) < 2i/ -|- 1, the Mellin transform of j^{x) is given by 
(38). Applying Parseval’s formula to these two Mellin transforms, we obtain, for 0 < 5 < 1, 


fOO 

/ jl{x)dx 

^0 

^ 02.-1 r (f) r + z.) r(i + ^)^ r (^) r(i + z.) 

2 «A-.oo r(i)r(i-f+ i.)r(i-f+ 2z.) r(i-^ + z.) 

2-^''+-^T(l + vf 1 rSG+^co 2-2‘r(t)r(i-t-f I/) T (§ -1) 

r(l) ^^A/2 -.oo T{i-t + n)r{i-t + 2,^)r{i- i + t + u) 


22iy+s — lj^]3 

r(^) 


G: 


3,3 


1 ,1 -f 1 -f 2iy 


The claim then follows, by analytic continuation, from Theorem 2.14 with n = 3 and k = 0. □ 

We note that, as in [BSW13], this Meijer G-function expression can be expressed as a sum of 
hypergeometric functions by Slater’s Theorem [Mar83, p. 57]. This is made explicit in (51). 

Equation (43) gives a hypergeometric expression for the even moments of a 3-step random walk. 
It was noticed in [BNSWll] that, in the case of planar walks, the real part of this hypergeometric 
expression provides an evaluation of the odd moments. These odd moments are much harder to 
obtain, and it was first proved in [BNSWll], based on this observation, that the average distance 
of a planar 3-step random walk is 


where 


W3(0; 1) = A+ 44 « 1-5746, 

TT^ A 

w^3(0;-i) = A^r‘^(i) 


(49) 

(50) 
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In the sequel, we generalize these results from two to arbitrary even dimensions. In particular, as 
explained in Example 3.12, we establish the transcendental nature of the odd moments of 3-step 
walks in all even dimensions by showing that they are all rational linear combinations of A and 
l/(7r^A). 

Theorem 3.10. (Hypergeometric form of H 3 at odd integers) Suppose that d is even, that 
is, V is an integer. For all odd integers s ^ —2v — 1, 


Wsii'^s) = Re 3 F 2 


/—s/ 2 ,—s /2 — v,iy + 1/2 
\ -f 1, 2u 1 



Proof. The case 1/ = 0 is proved in [BNSWll, Theorem 6]. We will prove the general case by 
induction on v. 

It is routine to verify that the hypergeometric function 


pz ^ T? -s ^,-s 


which, for even s, agrees with W 3 {i';s) by (43), satishes the contiguity relation 

(s -I- l)(s -I- 6z^ — l)F{i'; s — 1) -I- (!>v^F{v — 1; s -I- 1) — 2i/^F{v — 1; s -I- 3) = 0. 
On the other hand, it follows from Theorem 3.9 and Slater’s Theorem [Mar83, p. 57] that 


W3(.-s) = r(-.-^)r(. + i)^r(. + | + i) 

27rr(-f)r(2z/-b^) U-b^,2z/-b^ 


1 l-iy I 

2 ’ 2 ^’2 


2-+^-r(z. + i)r(zz + ^) ^ 

VttE (2z/+ f + 1) iy + l,-iy-^ 


(51) 


In that form, it is again a routine application of the holonomic systems approach [Zei90] to derive 
that 


0 = (s + l)(s + 6z/-l)IT3(zz;s-l) 

- 1; s -f 1) - 2 v^W 3 {v - 1; s -b 3). 


(52) 


Since this relation matches the relation satisfied by F[v, s), and hence Re F[v, s) when v and s are 
real, the general case follows inductively from the base case v = Q. □ 


Remark 3.11. The coefficients of the hypergeometric functions in (51) can be expressed as 

t{v-s) 22 '^-! 


Cl := 


/ 2i/+s+l \ 
\z>+(s+l)/2/ 


- (2-+i)e)e:/r)r,::r^ 


and C2 := (J+^/2)j respectively. Here, the factor 


T{y,s) = 


1 


cos(7rzz) cot(7rs/2) — sin(7ri/) 


is ±1 for half-integers v, and ±tan(^) for integers v. 
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Example 3.12. (Odd moments W 3 {i'; ■) in even dimensions) The planar case of Theorem 3.10 
was used in [BNSWll] to prove that the average distance to the origin after three random steps in 
the plane is given by (49). It is a consequence of (48), extended to complex k, that all planar odd 
moments are Q-linear combinations oi A = ^j^r®(l/3), dehned in (50), and 

The dimensional recursion (52) used in the proof of Theorem 3.10 shows that this observation 
extends to all even dimensions. For instance. 

Moreover, the average distance to the origin after three random steps in four dimensions is 

with similar evaluations in six or higher even dimensions. o 


Example 3.13. Theorem 3.10 does not hold in odd dimensions, in which the involved quantities 
can be evaluated in elementary terms. For instance, in the case of dimension 3, 


Wsi^is) 


1 3®+3 - 3 

4 (s + 2) (s + 3) 


while for integer s > 0 


3 


F2 


only agrees for even s. 



-s/2-1/2,1 

3/2,2 


1 3 «+ 3 - 2 -(- 1 )« 
4 (s + 2)(s + 3) 


o 


Example 3.14. (Eirst derivative of in even dimensions) By differentiating the 

hypergeometric representation of lF3(:^;s) in (51), it was shown in [BSWZ12, Examples 6.2 and 
6.6] that 

where the derivatives are with respect to s. It follows from differentiating (48), extended to complex 
k, that all derivatives 2k) lie in the Q-linear span of 1, ^ Cl(^) and Then, differentiating 

(52), we find that, indeed, for all integers > 0, the derivatives W/(z^; 2k) can likewise be expressed 
as 

W^{v]2.k) =ri +r2— +r3-Cl(^), 

TT TT O 

with rational numbers ri,r2,r3. 

Moreover, in the case of lF/(z/;0), a slightly more careful analysis reveals that r3 = 1. While 
we omit the details, we note that this can be seen, for instance, by evaluating 14^3(1; 0) and then 
deriving, in analogy with (52), a functional equation relating W^{v, s), s) and W3{i'+2-, s). 

In four and six dimensions, we obtain, for example, 


W'(1;0) 


1 

2 


11^3 

167r 


ici 

TT 



W'(2;0) 




A special motivation for considering these derivative values is that, in the case of two dimen¬ 
sions, lF/(0; 0) is the (logarithmic) Mahler measure of the multivariate polynomial 1 -I- xi -f 0:2; see 
[BSWZ12, Example 6.6] or [BS12, Section 4]. o 
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Example 3.15. (Second derivative of 113 ( 1 ^; ) in even dimensions) The second derivative 
0), interpreted there as a higher Mahler measure, is evaluated in [BS12, Theorem 4.4] in the 

form 

W^3'(0;0) = y + ^Ls3(^y), (53) 

where Ls„ denotes the nth log-sine integral 


Ls„((t) 




AO. 


For alternative expressions of the log-sine integral in (53) in terms of other polylogarithmic con¬ 
stants, we refer to [BS12]. We have not been able to obtain an equally natural log-sine evaluation 
of 

-I ^27r p ‘ 2 t ^ 

W"( 0 ; 2 fc) = — y J |l + e*‘" + e*®|2'=log2|l + e“ + e*V^dw, (54) 

for integers fc > 0. We may, however, derive 


W^'{0;2) = 3W^'{0;0) + ^ cosujReJ^ log^ (l - 2sin e*®) d6»dw. 
By the methods of [BS12, §4], we arrive at 

W^ 3 '( 0 ; 2 ) = 3TT3'(0;0) J ^ 12(1 sin^ coswdw. 

We may now integrate by parts and obtain 


W"(0;2)=3W"(0;0)-1 


2^3 

TT 




log s ds. 


(55) 


(56) 


(57) 


Moreover, 




log s ds 


00 


E 

n=0 


dn 

3" n^ 


(58) 


where a„ is given by [Slol5, A025565] of the OEIS and counts “the number of number of UDU-free 
paths of n — 1 upsteps (U) and n — 1 downsteps (D)” with recursion 


(n — l)a„ — 2(n — l)a„_i — 3(n — 3)a„_2 = 0. 

(Alternatively, a„ = M„_i Mk-iOn-k with Mk the Motzkin numbers given in A001006.) 

Solving for the generating function 03 ( 0 :) of a„ and considering a 3 (a;)da:^ we finally arrive 

at 

W"(0; 2) = 3 W"(0; 0) - 3 ^(log 3 - 1) - E -Cl (E. (59) 

As in Example 3.14, all second derivatives 2k) in even dimensions may then be expressed in 

terms of 1 T 3 '( 0 ; 0 ) and 1 F 3 '( 0 ; 2 ) as well as the constants in Example 3.14. For instance, 

W"( 1 ; 0 ) = -|W'"( 0 ; 0 ) -h ;^VF"( 0 ; 2 ) - E " E d f E , 

^ ^ ^ 8 ^ ^ ^ 24 ^ ^ ^ 8 487r Gtt V 3 / 

and so on. o 
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3.3 Moments of 4-step walks 

Lemma 3.16. The nonnegative even moments for a A-step walk in d dimensions are 

k fk\ r2j+2v\ r2{k—j)+2iy\ 

Wa (v 2k) = ' 3 ' ' j ' ^ ' 

’ ) 2-^ (k-j+v\ 

j=0 V j / V j / V k-j ) 

Proof. We apply Corollary 2.21 with ni = 2 and n 2 = 2, to obtain 


(60) 




i)er) 


j+l') 
3=0 V 3 I 


W2im,2j)W2iiy;2ik - j)). 


Using the evaluation of lU2(t^; 2k) given by (36) then yields (60). 


□ 


Example 3.17. (Generalised Domb numbers) The binomial sums in (60) generalize the Domb 
numbers, also known as the diamond lattice numbers [Slol5, A002895], 


W4(0;2fc)=^ 

3=0 


" 'kV ( 23 \( 2 {k-j) 


jj \j J\ fc-J / 


(61) 


for fc = 0,1, 2,..., which have played an important role in dimension 2. Their ordinary generating 
function, 


Y,W 40 -, 2 k)x^ = - 


1 


fc=0 


— 16a; 
1 

1 — 16a; 


1 1 


108x 


2.U1 


6’ 3 
1 


(16a; — 1)3 


108a; 


(16a; — 1)3 


(62) 


was determined in [Rog09]. The final equation follows from Clausen’s product formula. 
In four dimensions, the recursive relation (30) combined with (32) yields 


fc+i 


Wi{l-,2k)=J2Nik + l,j)C,Ck-j+2, 

i=i 


where Ck are the Catalan numbers, as in (32), and 


N{k +1,3 + 1) 


j + 1 \j 


k + 1 
j 


are the Narayana numbers, as in Example 3.17. After developing some further properties of the 
moments, we illustrate in Example 3.25 that the ordinary generating function for the even moments 
W 4 {v, 2 k) can be expressed in terms of hypergeometric functions whenever the dimension is even, o 

Corollary 3.18. (Hypergeometric form of W 4 at even integers) For k = 0,1,2,we have 


W4{iy; 2k) 


^2fe+2i/^ 



4^3 


/ —k, —k — V, —k — 2v, V + 1/2 
\ V + l,2v + 1,—k — V + 1/2 
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We note that this hypergeometric function is well-poised [DLMF, §16.4]. 

Example 3.19. Applying creative telescoping to the binomial sum (60), we derive that the mo¬ 
ments WA{v]2k) satisfy the recursion 

{k -\-2v l){k -I- -I- l)(fc 4:1/ \)W 4 ^{v\ 2k -\- 2) 

= ((fc + s) + 2*^) (20 {k+^Y +80 {k+^)i/ + 48v^ + s) Wa{v- 2k) 

-04k{k + v){k + 2v)WA{v,2k-2). (63) 


As in Example 3.8, Wi{v]s) is analytic, exponentially bounded for Re s ^ 0 and bounded on 
vertical lines. Hence, Carlson’s Theorem 3.7 applies to show that the recursion (63) extends to 
complex k. o 


The following result is the counterpart of Theorem 3.9 and extends [BSW13, Theorem 2.8]. 


Theorem 3.20. (Meijer G form of W 4 ) For all complex s with Re (s) > —4i/ — 2 and dimensions 
d'+2, 

4- V + Y ^ 1, — i/, 1 + v, 1 + 2i/ 


W4iv,s) = 


T I 




-g: 


4,4 

2,2 


-i/,i + i/,i 


Proof. The proof is obtained along the lines of the proof of Theorem 3.9. This time, Parseval’s 
formula is applied to the product jf,{x) = jl(x)j^{x). □ 

Remark 3.21. We note that, as in the case of (51) for three steps, the Meijer G-function in 
Theorem 3.20 can be expressed as a sum of hypergeometric functions, namely 


W4{t/;s) 


di ■ 4 F 3 


+d2 ■ 4 F 3 


Y h ~ 4~ 2^" -|- I -f 1 

^i/+ ^,2i/+ 

u + l,v + 2 ,- 1 / — 


(64) 


by Slater’s Theorem [Mar83, p. 57]. The coefficients di and d^ of the two hypergeometric functions 
are 


di 


2iu+s+iY{y + l)T{2i/ f -f 1) 

V^r(3i/ + ^) 


Cl, 


d2 ■— C 2 , 


where the coefficients Ci and C 2 are as in the three-step case, see Remark 3.11. 


Example 3.22. (Odd moments W 4 {i/; •) in even dimensions) It was shown in [BSW13, 
Section 2.3 & 3.1] that the average distance to the origin after four random steps in the plane, as 
well as all its odd moments, can be evaluated in terms of the elliptic integrals 


A 

B 


iri /5111111 

-A i^'(fc)2dyfc=-7F6 M’f’2’2-2-2>2 


[ k^K'{k)^dk 

Jo 


Stt 


7-^6 


7 3 3 1 1 1 1 
4’ 2’ 2’ 2’ 2’ 2’ 2 


3 

4’ 


2 , 2 , 2 , 2,1 
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We note that each of the -jFq hypergeometric functions may alternatively be expressed as the sum 
of two 6^5 hypergeometric functions. Then, in two dimensions, 

Wa{Q-,-1) =AA, W4(0;l) = 16A-48B, 

and it follows from (63), extended to complex values, that all odd moments are indeed rational 
linear combinations of A and B. In order to generalize this observation to higher dimensions, we 
claim that 


0 = 3(s + 2)(s + 6 i^)(s + 8 i^)(s + — 2 )W 4 (^; s) 

+256:/3(s + - 1; s + 2) - + 32iy - 6 )W 4 (:^ - 1; s + 4). (65) 


This counterpart of (52) can be proved using the holonomic systems approach [Zei90] applied to 
the hypergeometric form (64). We conclude that, in any even dimension, the odd moments lie in 
the Q-span of the constants A and B, which arose in the planar case. For instance, we hnd that 
the average distance after four random steps in four dimensions is 


W"4(l;l) 


3334144 ^ 11608064^ 

165375 ^ ~ 165375 ^ 


and so on. 


o 


Example 3.23. (First derivative of Wi(y] •) in even dimensions) In continuation of Exam¬ 
ple 3.14, we recall from [BSWZ12, Examples 6.2 and 6 . 6 ] that 


W((0;0) 


7C(3) 

27r2 ’ 


W^(0;2) 


14C(3) 12 , , 

9 o “T O. 

TT^ TT^ 


where, again, the derivatives are with respect to s. These evaluations may, for instance, be ob¬ 
tained from differentiating the hypergeometric expression (64). Proceeding as in Example 3.14, we 
differentiate both (63), extended to complex values, and (65), we conclude that, for all integers 
k,iy ^ 0, 

Wi{v]2k) = 

TT^ 27r2 

with rational numbers ri,r 2 ,r 3 . Again, we find that ra = 1 in the case k = 0. In four and six 
dimensions, we obtain, for example. 


Wi(l;0) 


3 53 7C(3) 

4 Qtt^ 27r2 ’ 


Wi(2;0) 


13 48467 

M ~ 141757r2 


7C(3) 

27r2 


In analogy with the case of three-step walks, the derivative 1 T 4 ( 0 ; 0 ) is particularly interesting 
because it is the Mahler measure of the multivariate polynomial 1 + xi + X 2 + x^. o 


Example 3.24. (Second derivative of W 4 ^{v\ •) in even dimensions) In [BSWZ12, (6.5), 
(6.10)] evaluations for the second derivatives 1 F 4 ( 0 ; 0 ) and 1 F 4 ( 0 ; 2 ) are given in terms of poly- 
logarithmic constants. It follows from these evaluations and the functional equation (63) that the 
derivative values 1 F 4 '( 0 ; 2k) all he in the Q-linear span of 


log 2 , —, 

TT^ 


1 log 2 C(3) 


log" 2 


C(3)log2 Li4(l/2) 
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where Li„(z) := is the polylogarithm of order n. Indeed, we realize from the dimen¬ 

sional recursion (65) that the same is true for the values W!^{v,2k) in all even dimensions. For 
example, 


n^ 253 , 5IF^(0;0) 2391^^(0; 2) 26IF^'(0;0) , 53W^'(0;2) 

+ 648 27 + 108 

25 1 2 , 2 o , 1193 1061og2 , 21 C(3) 

= -^-5^ -log + 

+ ^ + 24^MV2) . 

TT^ 7r2 7^2 

The number of (presumed) transcendental constants can be somewhat reduced when working in 
terms of Kummer-type poly logarithms, as illustrated in [BS12, Theorem 4.7]. o 

Example 3.25. (Ordinary generating function for even moments with four steps) In 

(62) we noted that the ordinary generating function of the moments 1 F 4 ( 0 ; 2k) has a concise hyper¬ 
geometric expression. It is natural to wonder if this result extends to higher dimensions. 

Combining the recursive relations (63) and (65), we are able to derive that the ordinary gener¬ 
ating function of the moments W 4 (z/; 2k), when complemented with an appropriate principal part 
(as in Theorem 3.6), can be obtained from the corresponding generating function of — 1; 2k) 
as well as its first two derivatives. Because the precise relationship is not too pleasant, we only 
record the simplified generating function, 


1 1 00 

^ ^ ^ W 4 (l; 2k)x'^ = (32a: - 7)^^ - (4 :e - 1) 


n—0 


that we obtain in dimension 4. Here, 


(32a: -f 3)FoFi - 16a:^ -f 10a: -f - F^ 


( 66 ) 


F. := 


2 F 1 


2 ■ 3'^x(16a: — 1)1“'^ dx^ 
which, by the differentiation formula [DLMF, (16.3.1)] 


1 1 

6’ 3 
1 


108x 


(16x — 1)3/ ’ 


^ F 


ai 


, ,u,p 


bi,...,bg 


(ai)n ■ ■ ■ (Op) 


(bl)n ‘ ' ' (bq)n 


PJn p 
pXq 


Ul -f 1, . . . , dp -f 1 

+ 1 , . . . , ^5 + 1 


(67) 


can be expressed as generalized hypergeometric functions. It would be nice if it was possible to 
make the general case as explicit as we did in Theorem 3.6 for three-step walks, but we have not 
succeeded in doing so. o 


3.4 Moments of 5-step walks 

As in the planar case, as well as in many related problems, it is much harder to obtain explicit 
results in the case of five or more steps. This is reflected, for instance, in the fact that an application 
of Corollary 2.21 with ni = 3 and 712 = 2, and appealing to (39), results in a double (and not single) 
sum of hypergeometric terms. 
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Lemma 3.26. The nonnegative even moments for a 5-step walk in d dimensions are 


W5{i^;2k) = 


k 

E 

j-o 

k 

E 

i=o 


}jr/) czr) 


-n’3W;2j) 


k\ (T) c'TT') y- ''A (’tne-fn 


i){‘T) (‘;'r) (■+■') 


(68) 


Example 3.27. As in the case of three and four steps, we can apply creative telescoping to the 
binomial sum (68) to derive a recursion for the moments Wf,{i']2k). In contrast to the three-term 
recursions (48) and (63), we now obtain a four-term recursion, namely 


[k -\-2v -\- l){k -I- -I- l){k ii/ l)(fc 5iy -k l)W^{y] 2k 2) (69) 

= a {v^k W^{v\ 2k) — b{v\ k)W^{v, 2k — 2) -\- 225k(k — l)(fc -\- v){k — 1 -|- 2k — 4), 

where 


a{v,m) := 35771“^-|-350j^m^-|-(1183 j^^-|-^) 

-b (I540z/^ + ^)iym+ (600i/^ -b -b , 
b{L'-k) := /fc(A:-b7/)(259fc2-bl295fc7^-b 14507^^-b 26). 

As in Example 3.8, s) is analytic and suitably bounded for Re s ^ 0, so that we may conclude 

from Carlson’s Theorem 3.7 that the recursion (69) extends to complex k. o 

Example 3.28. (Dimensional recursion for 145 ( 77 ; s)) Creative telescoping, applied to the 
binomial sum (68) for 144 ( 77 ; 2fc), allows us to derive the following more involved counterpart of the 
dimensional recursions (52), (65) for three and four steps. 

0 = 3(s + 2)(s + 4)(s-b 277 -b 2)(s-b 877 )(s- b IO 77 — 2)(s-b 1077 ) 144 ( 77 ; s) 

-45077^(5 + 4)(s + 277 -b 2 ) 144(77 - 1; s + 2) 

+ 477 ^ 0 ( 77 ; 5 ) 14 ^ 5(77 — 1; s -b 4) — 277'^6(77; s)144(77 — 1; s + 6), (70) 

where 

0 ( 77 ; s) ;= lOTs^ + 2(44577 -bl52)s-b 2(550772-b 116577 - 78), 

6 ( 77 ; s) := (s + 477 -b 2)(13s + IIO 77 — 16). 

This recursion is first obtained for nonnegative even integers s, and then extended to complex values 
using Carlson’s Theorem 3.7. o 


4 Densities of short random walks 

4.1 Densities of 2-step walks 

We find an explicit formula for the probability density ^ 2 ( 77 ; x) of the distance to the origin after two 
random steps in by computing the Bessel integral (6). An equivalent formula is given in [Wanl3, 
Corollary 4.2], which exploited the fact that the probability density is essentially the inverse Mellin 
transform of the moments which are evaluated in Theorem 3.2. 
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Lemma 4.1. The probability density function of the distance to the origin in d 2 dimensions 
after 2 steps is, for 0 < x < 2, 


P2{v;x) 





(71) 


Proof. It follows from (17) that 


P 2 {v,x) = v\2''x 


v+l 


P '' Ju{tx)Jf,{t)At. 


From [Wat41, Chapter 13.46], we have the integral evaluation 



t^ ''J,j{at)J,,{bt)Ju{ct)<Tt 


{abcYV{u + I/2)T {1/2)' 


(72) 


assuming that Re {v) > —1/2 and that a, b, c are the sides of a triangle of area A. In our case, 


A = f 



and therefore 


P2Y;x) 


which is equivalent to (71). 




r(:. + i/2)r(i/2)' 




2\<y-i/2 


□ 


Note that (71) reflects the general fact, discussed in Section 2.2, that the densities Pn(y\x) are 
piecewise polynomial in odd dimensions. 

For comparison with the case of 3 steps, we record that the density P 2 {v] x) satisfies the following 
functional equation: if F{x) := P 2 (v-, x)/x, then 


F{x) = F[y^-x^^. 


(73) 


Remark 4.2. The probability density in two dimensions, that is 


P2{Yx) 



is readily identihed as the distribution of 2| cos0| with 0 uniformly distributed on [0, 27r]. In other 
words, if X = {Xi,X 2 ) is uniformly distributed on the sphere of radius 2 in then p 2 { 0 ;x) 
describes the distribution of jXij. It is then natural to wonder whether this observation extends to 
higher dimensions. 

Indeed, if X = (Xi, X 2 ) is uniformly distributed on the sphere of radius 2 in where X is 

partitioned into Xi,X 2 G then p 2 (i';x) describes the distribution of jXij. Details are left to 

the interested reader. In this stochastic interpretation, the invariance of (73) under x 1 —>■ Vd — 
is a reflection of the fact that jXij and IX 2 I = -^/d — jXiP share the same distribution. 
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4.2 Densities of 3-step walks 

It was shown in [BSWZ12] that the density psiO; x) of the distance to the origin after three random 
steps in the plane has the closed form 


P3(0;a;) 


^ 1 (3 + cc2)3 ) ’ 


(74) 


valid on the interval (0,3). We next generalize this hypergeometric expression to arbitrary dimen¬ 
sions. In order to do so, we need to establish the behaviour of x) at the end points. 

To begin with, we use information on the pole structure of the moments s) to deduce the 

asymptotic behaviour of P 3 {v;x) as cc —)■ 0+. 

Proposition 4.3. For positive half-integer v and x —>■ O’*', 

(75) 

V 37r (J 

Proof. As shown in Corollary 2.15, the moments s) are analytic for Re s > —d and the first 

pole has residue 

POO n oiy 

ReSs=-dW3{iy,s) = v\2'' / x^~''J^{x)dx=—^j^, 

Jo v37r(^) 

where the last equality is another special case of (72). The asymptotic behaviour (75) then follows 
because W 3 (^; s — 1) is the Mellin transform of psii^; x). □ 

On the other hand, to obtain the asymptotic behaviour of psii^; x) as x —>■ 3“, we have to work 
a bit harder than in the case of the behaviour (75) as a; —>■ 0“*". 

Proposition 4.4. For positive half-integer v as x ^ , 

./a p2i^oi^ 

P3{<^;x)^^jj;y{3-xrr (76) 

Proof. Using Theorem 2.9 together with the fact that, for x G [0,2], 


we find 




(2x)2'^+i 



(4- (l-k2Ax-kx^))'^-^/^ 
(1 -I- 2Ax -I- x2)3/2 


(1 - 


Observe that the upper bound of integration is 1 if x G [0,1], and (3 — x^)/(2x) for x G [1,3]. On 
X G [1,3], after substituting t = Vl + 2Ax -I- x^/x, we therefore find 


Po{^]x) 






(x + l)2)(tV 


(x-l)2)}"-i/2dt. 
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Note that the polynomial in the integrand factors into (2 — tx){tx — (x — 1)) times a factor which 
approaches 192 as x —>■ 3“ and t ^ 2/3. Hence, as x —>■ 3“, 

^ {( 2 -te)(te-(x-i))r-v 2 dt. 

(J ^ 

We now relate this integral to the incomplete beta function to find 

£j( 2 - fa)(fx - (i - l)))->/2dt = ■ 

Putting these together, we conclude that (76) holds. □ 

We are now in a position to generalize (74) to higher dimensions. 

Theorem 4.5. (Hypergeometric form for p^) For any half-integer ^ ^ 0 and x G [0,3], we 
have 

P3{iy;x) 2V3 3-3^ x^£9 - x")"^ /i,| 

X 3 + x2 ^ ^ 1 + ^ 

Proof. We observe that both sides of (77) satisfy the differential equation A 3 ■ y{x) = 0, where the 
differential operator A 3 is given by 

x(x^ - 9)(x2 - l)Dl - {2v - l)(5x^ - 30x^ + 9)i:)„ + 4x(x2 - 3){3v - 1)(2^ - 1), 

where Note that, in the case of P 3 (z^;x), this is a consequence of the functional equation 

resulting from (48). The indices of this differential equation at x = 0 are 0 and 2v. It follows from 
Proposition 4.3 that P 3 (v]x)lx is the unique solution j/(x), on (0,1), such that 

v37r (J 

as X —>■ O'*'. Since this property is satisfied by the right-hand side of (77) as well, it follows that 
(77) holds for all x S (0,1). 

Similarly, to show (77) for all x G (1, 3), we use the fact that the differential equation has indices 
0 and 2 i 2 at X = 3 as well. In light of Proposition 4.4, p 3 {v]x)/x is the unique solution y{x), on 
(1,3), such that, as x —>■ 3“, 

1 (A2voy 

Again, it is routine to verify that this property is also satisfied by the right-hand side of (77). 
Hence, it follows that (77) holds for all x G (1, 3). □ 

As a consequence, we have the following functional equation for the probability density function 
P 3 {i'] x). The role of the involution x 1 —>■ \/4 — x^ for 2 steps, see (73), is now played by the involution 


x"^(9 — x^y 


(3 + x2)2 


(77) 


32 











Corollary 4.6. (Functional eqnation for p^) For any half-integer i/ 0 and x S [0,3], the 
function F(x) := p'i(y\x)/x satisfies the functional equation 


F{x) = 


l-\- X 


Gv-2 


F 


3 — X 

1 + a; 


(78) 


Proof. The hypergeometric right-hand side of (77) clearly satisfies the functional equation (78). □ 


It would be very interesting to have a probabilistic interpretation of this functional equation 
satisfied by the densities p^{i>]x). 

Remark 4.7. We note the following relation between the functional equations for the two- and 
three-step case in (73) and (78), respectively. Namely, if 


_ x + y 
2 


+ Ij 


x-y 
2i ’ 


then the relation y = translates into = 4 — X^. It is natural to wonder whether this 
observation might help explain the functional equation (78) for three-step densities. 


Example 4.8. As a consequence of the contiguity relations satisfied by hypergeometric functions, 
we derive from (77) that the densities p 3 (i';x) satisfy the dimensional recursion 


P3{iy1; x) 


&{. 2 y+i)%tl\)(Su+ 2 ) - 3)(a;^ - 6a; - 3){x^ + 6x - ?,)pz{v, x) 

12(2j/-i)(2i/+i)(3!/-i-i)(3i/-i- 2) ^ ~ ~ Ps~ 1 ! a;), 


where v ^ 1. 


o 


Example 4.9. It follows from the hypergeometric formula (77) that the densities pz{v\x), for 
V > 0, take the special values 

^ 47r2 12 (2z/)!(3^)!‘ 

In particular, P 3 (^; 1) G Q in odd dimensions, and psiiy; 1) G Q • in even dimensions. o 

Example 4.10. From (19) and the functional equation of Corollary 4.6, for > 1, we learn that 

1) = 1) - - l)(6z/^ -\-5v-\- 2 )p3{i2; 1), 

but learn nothing about 1 ). o 


4.3 Densities of 4-step walks 


It is shown in [BSWZ12] that, in the planar case, the probability density of the distance to the 
origin after four steps admits the hypergeometric closed form 


P4(0;a:) = 


2 Vl6 - a;2 


Re 3 F 2 


111 
2 ’ 2 ’ 2 
5 Z 
6 ’ 6 


(16 — X^) 


2\3 


108a:4 


(79) 


In this section we obtain a higher-dimensional analog of this evaluation by demonstrating that 
hypergeometric formulae can be given for the densities P 4 {i'',x) in all even dimensions. 
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Example 4.11. Since the Mellin transform of the density p 4 ,{v\x) is given by the corresponding 
probability moments Wi{v]s — 1), the recursion (63) for the moments translates into a differen¬ 
tial equation for the density. We refer to [BSWZ12] for details. We thus find that pi{v]x)/x is 
annihilated by the differential operator 

x^{x^ - 16)(x2 - 4)Dl - 3(2iy - l)x‘^{3x^ - AOx^ + 64)DI 
+(2iy - l)x(i52i^ - I9)x^ - (416i/ - 152)x^ + 256iy - 64)D^ 

-8(4i/ - 1)(3 j^ - l)(2z^ - l)x‘^{x‘^ - 4), (80) 

where We conclude that, in the planar case, p4(0;a;) satisfies the differential equation 

0 = a;^(a:^ — 16)(a;^ — 4)p4"(0; x)-I- 6x^(x'^ — 10a;^)p4(0; x) 

-|-x(7x^ — 32x^ -I- 64)p4(0; x) + (x^ — 64)p4(0; x), (81) 

which agrees (up to a typo there) with [BSWZ12, (2.7)]. 

The differential equation (81) is the symmetric square of a second order differential equation and, 
moreover, admits modular parametrization [BSWZ12, Remark 4.11]. These ingredients ultimately 
lead to the hypergeometric closed form (79). The differential equation associated with (80), on the 
other hand, is a symmetric square of a second order differential equation only in the cases ly = 0 
and iz = 1/2. o 

Example 4.12. In three dimensions, that is when zz = 1/2, the density is 

P4(1/2;x) = :^(8 - x2 + 2(x - 2)lx - 2]), (82) 

16 

for X G [0,4]. This is equivalent to (24) and may be derived directly from (6) or from Theorem 2.11. 

o 


Example 4.13. Basic Mellin calculus connects the asymptotic behaviour of Pn{v]x) as x —>■ 0+ 
with the nature of the poles of Wn{v]s) in the left half-plane. For instance, from the explicit 
information in Example 2.17 on the poles of W 4 (l; s) at s = —4 and s = — 6 , we conclude that 


P4(l;a;) = -^x^ G- 
TT^ 


IGtt^ 


91og(2) \ 


+ ^ log(a^)x'' + O(x^) 


as X —>■ 0 +. 


o 


The following result connects the 4-step density in d dimensions with the corresponding density 
in d — 2 dimensions. In particular, using (79) and (82) as base cases, this provides a way to obtain 
explicit formulas for the densities p 4 (iz;x) in all dimensions. 

Theorem 4.14. (Dimensional recursion for P 4 ) For 0 ^ x ^ 4 and any half-integer ly ^ 0, 


3{2iy + l)(3zz -b l)(3z^ -b 2)(4z/ -b l)(4i/ -b 3) 


P4(j^ + l;a;) 


64(z/ -b 1)^ 

= —a{iy; x^/ 8 )p 4 (z/; x) -b xb{v\ xP' jFdfp^fy^ x) — c{v] x^/‘F)pf{v, x), 


(83) 
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where 


a{v]x) := (4i/-l)(6z/-l)x^ + 2(100z/^-23z/-l)x3 + 2(2i/ + 3)(12z/ + l)a;2 

-2{60v^ + Uv + l)x + {2v + l){Av + 1), 
b{v]x) := (lOz^ — 3)a:"‘+ 5(12:/— l)a;^ — — l)x^ — 20:/a: + 6z/+ 1, 

c{v;x) := 4a;(x — 2)(2a; — l)(a;^ + 5a; + 1). 

Proof. We have already observed in (63) and (65) that the moments W 4 (z/;s) satisfy a functional 
equation connecting W 4 (:/; s), W 4 (:/; s + 2), W 4 (:/; s + 4), as well as a functional equation relating 
W 4 {iy + 1; s), W 4 (:/; s + 2), W 4 {i'; s + 4). The usual Mellin calculus translates the first of these two 
into the third-order differential equation recorded in (80), while the second is translated into a more 
complicated equation involving seven terms ranging from P 4 (z/; a;) to {v + \]x). With assistance 
of the Mathematica package HolonomicFunctions, which accompanies Koutschan’s thesis [Kou09], 
we compute a Grobner basis for the ideal that these two relations generate and use it to find the 
claimed relation involving P 4 {v + 1; x) as well as P 4 {v\ x) and its hrst two derivatives. □ 


Example 4.15. (Hypergeometric form for p 4 (l;a;)) By combining (79) with Theorem 4.14 
and (67), we conclude, for instance, that the 4-step density in four dimensions, for x G (2,4), can 
be hypergeometrically represented as 


P4{P,x) 


(16-a;2)5/2 

{2A'kYx 


-{x^ + SfGQ + 


r{x‘^) 
7!! x'^ 


Gi + 


(16 — x^)^s(a;^) 
13!!(2/3)4x8 


(84) 


with (2n -I- 1)!! = (2n -|- l)(2n — 1) • • • 3 • 1 and 


Gx := 3 F 2 


^-l-A, ^-|-A, |-|-A 
i + A, g -b A 


(16 — a;^)^\ 
108a;4 j 


as well as 


r{x) := x^ + 55a;^ -b 1456a;3 -b 25664^2 - 90112x - 262144, 
s{x) := (x — 4)(a; -b 32)^(a;^ -b 40a; -b 64). 

As in (79), taking the real part of the hypergeometric functions provides a formula for p 4 (l;a;) 
which is valid for x € (0, 2) as well. The above hypergeometric formula also provides a modular 
parametrization (in a suitably generalized sense) for p 4 (l;a;). This is a consequence of 

P 4 ( 0 ; 8 ia;(r)) = /(t), 

TT 

which is proved in [BSWZ12, (4.16)] and which involves the modular function 

^ 77(2r)37?(6r)3 
ry(T)3?y(3r)3 

and the weight 2 modular form 

/(r) = ? 7 (r)? 7 ( 2 r)? 7 ( 3 T)? 7 ( 6 r). 

Differentiating this modular parametrization of P 4 ( 0 ; x), we find that P 4 (l; 8ix(T)) can be expressed 
in terms of modular quantities such as /'(r)/a;'(r). Then inductively we obtain the like result for 
higher even dimensions. o 
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Example 4.16. 

proves that 


(P4 (i/;2 ) in even dimensions) Motivated by [BSWZ12, Corollary 4.8], which 


P4(0;2) 



V3 


1^3(0;-1), 


or 7r/V3P4(0;2) = 1 ^ 3 ( 0 ;-!), we discover that, for integers i' 0, Tr/\^pi{v;2) is a rational 
combination of the moments 14^3(0; —1) and 14^3(0; 1). For instance. 


^P 4 (l; 2 ) = -y 1F3(0; - 1 ) + y VF3(0; 1 ) 


and 


73 


P4(2;2) 


6656 

715" 


W"3(0;-l) 


704 

"67 


4F3(0;1). 


To deduce the first of these from equation (83) takes a little care as p 4 (l;a;) is not differentiable 
at 2 and one must take the limit from the left in (84); likewise for higher dimensions (note that 
c{v; 2) =0). Theorem 4.14 specializes to 


(2z/ + 1)(3^' + 1)(3^' + 2)(4r' + 1)(4^' + 3)p4(iz + 1; 2) 

= 4(1/ + 1 ) 7 ( 61 / + l)(12i/ - 7)p4{i'; 2) - 3Q{6iy - l)p'4(i/; 2)]. 

On the other hand, specializing a corresponding relation among P 4 {v; x), P 4 {i' + 1; x), P 4 {v + 2; x) 
and P 4 {i'; x), which one obtains as in the proof of Theorem 4.14, we find 

(2iy - l)(2iy + l)(3iy - 2)(3iy - l)(3iy + l)(3iy + 2)(4iy + l)(4iy + 3)p4(iy + 1; 2) 

= -16(iy + l)^(2iy - l)(3iy - 2)(3iy - 1)(547 + l)p4(iy; 2) 

+576iy^(iy + 1)^(6:/ — 5)(6j/ — l)p4(i/ — 1; 2). 


o 


4.4 Densities of 5-step walks 

The 5-step densities for dimensions up to 9 are plotted in Figure 4. 

A peculiar feature of the planar density is its striking (approximate) linearity on the initial 
interval [0,1]. This phenomenon was already observed by Pearson [Pea06], who commented on 
P5(0;x)/x, between x = 0 and a; = 1, as follows: 

“the graphical construction, however carefully reinvestigated, did not permit of our 
considering the curve to be anything but a straight line... Even if it is not absolutely 
true, it exemplifies the extraordinary power of such integrals of J products [that is, (6)] 
to give extremely close approximations to such simple forms as horizontal lines.” 

Pearson’s comment was revisited by Fettis [Fet63], who rigorously established the nonlinearity. 
In [BSWZ12, Theorem 5.2], it is shown that the density satisfies a certain fourth-order differential 
equation, and that, for small x > 0, 


P5(0;x) = 

A;=0 


(85) 
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Figure 4: x) for ^ = 0, i, 1,..., | 


where 


(15(2fc + 2)(2/c + 4))2r5.fc+2 
= (259(2fc + 2)'‘ + 104(2fc + 2)2)r5,fc+i 

-(35(2fc + if + 42(2fc + 1)2 + 3)r5,fc + {2kfr^^k-i, ( 86 ) 

with initial conditions ^5 _i = 0 and 

f5,k = ReSs=_ 2 fc -2 W^5(0; s)- 


Numerically, we thus find that 

P 5 ( 0 ; x) = 0.329934 x + 0.00661673 + 0.000262333 x^ + 0{xf, 


which reflects the approximate linearity of P 5 ( 0 ;a;) for small x. 

By (18), the residue rs^o = P 5 ( 0 ; 0 ) can be expressed as P 4 ( 0 ; 1). The modularity of P 4 in the 
planar case, combined with the Chowla-Selberg formula [SC67], then permits us to obtain the 
explicit formula [BSWZ12, Theorem 5.1] 


V5 r(i)r(^)r(^)r(A) 

"5.0 - ^ ^4 


(87) 


Moreover, high-precision numerical calculations lead to the conjectural evaluation [BSWZ12, (5.3)] 


’’5,1 = 


13 


r5,o - 


2 1 


( 88 ) 


and the recursion ( 86 ) then implies that all the coefficients in (85) can be expressed in terms 
of rs^o- 

Theorem 4.17. (Residues of W 5 ( 0 ;s)) The conjectural relation ( 88 ) is true. 
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Proof. As noted after Corollary 2.15, we have 


ReSs =_4 W^il; s) = pi{l] 1 ). 

On the other hand, applying the dimensional recursion (70) for s) with v = \ and using the 

values W 5 ( 0 ; 0) = 1 and W^(f)\ 2) = 5, we obtain 

7 25 35 7 

ReSs=_4 s) = g “ ~ 

since Vb 5 ( 0 ;s) has simple poles only [BSW13, Example 2.5]. As noted in [BSWZ12, (6.2)], it is a 
consequence of the functional equation (69) that 

225r5,i = 26r5^o - 16 - 201^5(0; 0) + 41E^(0; 2). 


Combining the last three equations, we arrive at 


107 525 


Equation ( 88 ) is therefore equivalent to 


1 105 1 

134 1 ; 1 = -• 

6 167r4 ra^o 


(89) 


The equality (89) can now be deduced (at least in principle) from the hypergeometric formula for 
P 4 (l;a;), made explicit in Example 4.15, the modular parametrization of p 4 ( 0 ;a:) as well as the 
Chowla-Selberg formula [SC67]. □ 

In conclusion, we know that P 5 ( 0 ;a;) has a Taylor expansion (85) at a: = 0, which converges 
and gives its values in the interval [0,1]. Moreover, we have a recursive description of the Taylor 
coefficients and know that they are all Q-linear combinations of in (87) and l/( 7 r"*^r 5 _o)- All of 
these statements carry over to the 5-step densities x) in any even dimension. Since the details 
are unwieldy, we only sketch why this is so. 

Recall that the moments W^{v]s) satisfy the recursive relations (69) and (70). Indeed, there is 
a third relation which connects IT 5 (j^; s), 1 ^ 5 ( 1 ^; s -I- 2), W5{i'+1-, s), W^{v + 1-, s-|-2). As in the proof 
of Theorem 4.14, the Mellin transform translates these three recursive relations into (complicated) 
differential relations for the densities P 5 {v\ x). Assisted, once more, by Koutschan’s package Holo- 
nomicFunctions [Kou09], we compute a Grobner basis for the ideal that these three differential 
relations generate. From there, we find that there exists, in analogy with Theorem 4.14, a relation 

x^Pz{v -k 1; x) = Ap 5 (i/; x) + Bp'^{v, x) + Cp'f{v; x) + Dp'f'{v, x), 


with A, B, C, D polynomials of degrees 12,13,14,15 in x (with coefficients that are rational functions 
in v). 

We therefore conclude inductively that, for integers v, the density ^ 5 ( 1 ^; x) has a Taylor expansion 
(85) at a; = 0 whose Taylor coefficients lie in the Q-span of in (87) and l/( 7 r‘^r 5 _o). It remains 
an open challenge, including in the planar case, to obtain a more explicit description of P 5 (z/;x). 
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5 Conclusion 


We have shown that quite delicate results are possible for densities and moments of walks in ar¬ 
bitrary dimensions, especially for two, three and four steps. We find it interesting that induction 
between dimensions provided methods to show Theorem 4.17, a result in the plane that we previ¬ 
ously could not establish [BSWZ12]. We also should emphasize the crucial role played by computer 
experimentation and by computer algebra. One stumbling block is that currently Mathematical 
and to a lesser degree Maple, struggle with computing various of the Bessel integrals to more than 
a few digits — thus requiring considerable extra computational effort or ingenuity. 

We leave some open questions: 

• The even moments W„(0; 2k) associated to a random walk in two dimensions have combina¬ 
torial significance. They count abelian squares [RS09] of length 2k over an alphabet with n 
letters (i.e., strings xx' of length 2k from an alphabet with n letters such that x' is a permu¬ 
tation of x). As observed in Example 2.22, the even moments W„(l; 2k) are positive integers 
as well and we have expressed them in terms of powers of the Narayana triangular matrix, 
whose entries count certain lattice paths. Does that give rise to an interpretation of the even 
four-dimensional moments themselves counting similar combinatorially interesting objects? 

• As discussed in Example 3.17, in the case v = 0, the moments W4{i';2k) are the Domb 
numbers, for which a clean hypergeometric generating function is known. Referring to Ex¬ 
ample 3.25, we wonder if it is possible to give a compact explicit hypergeometric expression 
for the generating function of the even moments W 4 {i'; 2 k), valid for all even dimensions, as 
we did in Theorem 3.6 for three-step walks. 

• Verrill has exhibited [Ver04] an explicit recursion in k of the even moments W„(0;2fc) in the 
plane. Combined with a result of Djakov and Mityagin [DM04], proved more directly and 
combinatorially by Zagier [BSWZ12, Appendix A], these recursions yielded insight into the 
general structure of the densities Pn(0; a;). Eor instance, as shown in [BSWZ12, Theorem 2.4], 
it follows that these densities are real analytic except at 0 and the positive integers n,n — 
2, n — 4,... It would be interesting to obtain similar results for any dimension. 

• By exhibiting recursions relating different dimensions, we have shown that the odd moments 
of the distances after three and four random steps in any dimension can all be expressed in 
terms of the constants arising in the planar case. Is it possible to evaluate these odd moments 
in a closed (hypergeometric) form which reflects this observation? 

• In the plane, various other fragmentary modular results are (conjecturally) known for five 
and six step walks, see [BSWZ12, (6.11), (6.12)] for representations of IE5(0; 0) and IEe(0; 0), 
conjectured by Rodriguez-Villegas [BLRVD03], as well as a discussion of their relation to 
Mahler measures. Are more comprehensive results possible? 
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